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We study the ground state properties of a one-dimensional Bose gas with iV-body attractive 
contact interactions. By using the explicit form of the bright soliton solution of a generalized 
nonlinear Schrodinger equation, we compute the chemical potential and the ground state energy. 
For N = 3, a localized soliton wave-function exists only for a critical value of the interaction 
strength: in this case the ground state has an infinite degeneracy that can be parameterized by the 
chemical potential. The stabilization of the bright soliton solution by an external harmonic trap is 
also discussed, and a comparison with the effect of TV-body attractive contact interactions in higher 
d , dimensions is presented. 

^t- : 

C\J . I. INTRODUCTION 

One-dimensional quantum systems have always attracted a lot of theoretical interest: their exact solutions give, in 
fact, useful insights on the role of the interactions and other non-perturbative features JlL 3, [3 , H . In this respect, 
the growing ability to realize and manipulate one-dimensional Bose gases @, I, 0, S I, M, Mill O has provided 
^ ' an highly controllable experimental counterpart to these theoretical achievements. 

Quasi-one-dimensional Bose gases are obtained by using a cigar-shaped external trapping potential, elongated in 
a direction, with the other degrees of freedom frozen due to the presence of a tight transverse confinement. In the 
experiments, several variants of the interacting Bose gas in one dimension can be implemented: an optical lattice 
can be added to detect the Mott-superfluid transition in one dimension I14L the effective one-dimensional interaction 
can be tuned [l5[ to observe a Tonks-Girardeau gas of ultracold atoms [al§] , or the effect of the temperature can be 
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studied [13j. Important tools that permit to further control the pro per ties of low-dimensional Bose systems are the 
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tuning of an external magnetic field near a Feshbach resonance UflUa, and in perspective, the implementation of 



tuning oi an external magnetic neia near a fesnbacn resonance (la . uyij . ana, m persp 
the recently proposed schemes to engineer effective three-body interactions [18|, [l9|, [2(| ■ 

The technique of Feshbach resonances permits to change the sign of the scattering length: by switching from 
repulsion to attraction, i.e. from positive to negative scattering length, the homogeneous 1-D Gross-Pitaevskii equation 
(GPE) admits a solution corresponding to a localized wave- function, the so-called bright soliton 0, [T 5 , [l 7l . Bright 
matter- wave solitons were created both in Bose- Einstein condensates of 7 Li p, [Tfl and 85 Rb atoms 1211 . Various 
localized states has been also produced in quasi-one-dimensional geometries (for reviews see the book |22j|). 

With attractive two-body interactions, a crucial role is played both by the dimension of the system and the trapping 
potential. In three dimensions, for instance, homogeneous attractive bosons are unstable against the collapse, but the 
presence of an external harmonic trap can stabilize them: the critical value of the interaction coupling that gives rise 
to the collapse can be obtained from the GPE [l7| , and the critical particle number is given by ~ Mr \ a | / a osc 
t-h ■ where N't is the total number of particles, a < is the scattering length and a osc is the harmonic oscillator length 
[23L I24I I25I l2q] . In the one-dimensional case, the bright soliton solution is the ground state of the homogeneous GPE 
with negative scattering length. Furthermore, the GPE ground state energy is in agreement, in the thermodynamic 
limit, with the ground state energy obtained by Bethe ansatz for the attractive one-dimensional Bose gas [27| (see 
more in Section II). 

In this paper, motivated by the recent papers [Hi , |20| in which different schemes have been proposed to realize 
effective tunable three-body interactions, we consider an attractive three-body contact potential and, more generally, 
a TV-body contact interaction. We consider the limit of large number of particles, Mr >> 1, with the constraint 
cAf^ N 1 ' = const (c being the strength of the A-body interactions) so that the energy per particle is finite. Since no 
Bethe solution is available in the general case of A-body interaction, we employ an Hartree approximation to study 
the problem in the limit mentioned above . This means that the ground state energy is estimated by using the bright 
soliton solution of a generalized mean-field GPE equation. As we will show, the A = 3 is a special case: for this 
value, in fact, a localized soliton wavefunction exists only for a critical value of the interaction strength and has an 
infinite degeneracy. The stabilization of this bound state can be cured by putting the system in an external harmonic 
trap. The variational approach, that we will also employ, reveals the tendency of the higher body interactions to 
become more unstable in higher dimensions. It is worth stressing that the case we are considering does not consist of 
a A-body interaction added to the 2-body interaction of the Bose gas: we are interested, in fact, to the effect of the 
A-body in its own, since the coefficient of the two-body interaction can be tuned to be zero [20| . 

The plan of the paper is the following: in Section II we introduce the Hamiltonian corresponding to A-body contact 
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attractive interactions and we write the (mean-field) generalized GPE. The familiar case N = 2 is briefly recalled. 
In Section III the bright soliton solution for the homogeneous limit is obtained by using a mechanical analogy with 
a fictitious particle moving in a potential, and its properties are investigated. The comparison with the numerical 
results confirms that for N < 3 this is the ground state of the generalized GPE, as expected. The ground state energy 
by varying N is also determined. In Section IV we consider the effect of an harmonic trap: using a variational ansatz 
for the ground state we determine the critical value of the interaction needed to stabilize the bound state. In Section 
V there are our conclusions. 



II. TV-BODY ATTRACTIVE CONTACT INTERACTIONS 

The general quantum Hamiltonian for an homogeneous one-dimensional Bose gas with TV-body interactions 
V(x\,- ■ ■ ,x N ) is 



( Xl ) ■ ■ ■ W{x N )V( Xl , ■ ■ • , x N )V(x N ) ■ ■ ■ V( Xl ) , (1) 



where ^(x) is the bosonic field operator. The Lieb-Liniger Hamiltonian for the interacting one-dimensional Bose 
gas [28[ has the kinetic term plus a density-density term involving pairs of particles interacting via a contact two- 
body potential; this corresponds to TV = 2 and V{x\,X2) — Vq5{x\ — X2)' Vq positive (negative) corresponds to 
repulsion (attraction) between the bosons. The low-energy properties of the Lieb-Liniger model can be studied by 
the Luttinger liquid effective description [29[ obtained by bosonization [3(j (a general discussion of the correlation 
functions is presented in [3l1|). 

For TV-body attractive contact interactions we set V(xi, • • • , x^) = — cj^^ 1 5{xi — Xi + i) (c > 0): the Hamiltonian 
(P) reads then 



H = j dx¥{x) (~^^) ^( x ) dx ^( x ) ■ ■ ■ * t (-^)*(a;) • • • * 

In the Heisenberg representation, the equation of motion for the field operator is given by 



(x) . (2) 



iH -dt 



h 2 d 2 ~ / - x\ M-l / „ \ AT-1 A 



For N = 2, the corresponding Lieb-Liniger model is integrable and the ground state energy E can be determined by 
Bethe ansatz (32j: the final result is given by 

E_ mc 2 (A/g - 1) 

where N't is the total number of particles. For large Mr, from ((4]) it follows that one has to keep the product 
cMr — const in order to have a finite ground state energy per particle. Using the integrability of the N = 2 model, 
the correlation functions of the attractive one-dimensional Bose gas at zero temperature were recently calculated in 

m. 

For the three-body problem (N = 3), no Bethe ansatz solution is available, except for a more complicate double-<5 
function Bose gas which can be mapped in a one-dimensional anyon gas [3J|. Hence, to estimate the ground state 
energy E we propose here to employ a mean-field (Hartree) approach: in this approach, the ground state energy is 
given in terms of the ground state energy of a generalized GPE. The same procedure will be employed for other values 
of AT. 

Before we start the discussion of the general TV-body case, let's briefly remind how this task can be successfully 
done for N = 2 (27j . First of all, in the mean-field approximation the ground-state wavefunction is written as 

N't 

i>Gs{xi, ■■■ ,XU T ) OC JJVo(^i) , (5) 

i=l 

where the function ipo(x) is the ground state of the time-independent homogeneous GPE, i.e. the nonlinear Schrodinger 
equation (NLSE), given by 

h 2 d 2 
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where [i is the chemical potential and the normalization is given by J dx | ipo \ 2 — Nt- The energy is expressed as 

h 2 d 2 c 



E G p = J dxTpoix) 
The static bright soliton solution of © is given by 



M*) ■ ( 7 ) 



- ^c^b - (8) 



with k — ijicNt /2h 2 and N = (1 / '2) W 'mcNr / 'h 2 . Substituting this expression in ([7]) one gets 

E GP _ mc 2 N% 
~Nt ~ 24ft 2 ' 



(9) 



i.e., the exact result |(U) apart terms esc \fN T . A comment is in order: in the homogeneous one-dimensional interacting 
case there is, strictly speaking, no condensate. However the condition cNt — const implies that, for large Nt, the 
coupling constant should scale to zero, c — > 0: hence, we are in a weak coupling regime where the mean-field GPE 
is expected to give reasonable results. In a similar way, for c < (repulsive interaction) the comparison between 
the exact and the GPE ground state energy shows that the latter gives the correct behaviour for c — > while the 
Bogoliubov approximation gives the exact first-order corrections for small | c | [28| . 

Based on the analysis above, for general N and in the limit c^Owe expect that a reasonable description of both 
the ground state properties and the low-energy dynamics is given by the mean-field generalized homogeneous GPE 

ih^§^> = (-^^--c\^{x,t)A^{x,t) , (10) 



dt \ 2m dx 2 

where the nonlinearity degree a of the generalized nonlinear Schrodinger (GNLSE) equation is related to N by 

N=~ + l. (11) 
The mean-field ground state is given by the time-independent GNLSE equation 

"|^^2 - <#o(z)| a ) M*) = , (12) 

where, as before, /i is chemical potential and ipo is normalized to the total number of particles N't, i.e. J dx\ipo(x)\ 2 = 
N T - 

With N > 2 integer, a is an even integer; however, in eqn. (|12[) a can take any real positive value and, in the 
following, we will consider this general case. In this respect let's comment that the axial dynamics of a Bose-Einstein 
condensate induced by an external potential with cylindrical symmetry in the transverse directions can be studied 
by introducing an effective one-dimensional GPE equation with a = 1 [35| and that for Bose-Einstein condensates in 
one-dimensional optical lattices the effective equation has a value of a that depends on the details of the trapping 
potentials and it is, in general, a non-integer value [36| . 

III. GROUND STATE OF THE GENERALIZED NONLINEAR SCHRODINGER EQUATION 

In the following we will study the attractive iV-body problem in the thermodynamic limit, defined by Nt —> oo, with 
the product G = cN^ 2 kept fixed. This will ensure the energy per particle of the GNLSE bright soliton to be finite. 
In dimensionlcss units, rescaling the wave function ipo — > y/Nr^o, cqn. (p~2|) reads 

1 d 2 \ 

2Q-2 -amwr I vfo(x) = fr*fo(x) , (13) 

where g and jl are the dimensionless versions of G and fj, respectively. For the ground state of this equation we look 
for a real solution, with the normalization condition 

ip%(x)dx = 1 . (14) 



FIG. 1: Typical shape of the potential V(tpo) f° r negative values of fi. ±A are the inversion points of the motion. 



Obviously, once a static solution ipoi x ) of ec l n - (|13|) has been found, the corresponding soliton wave solution with 
velocity v is given by 

$o(x,t) = - vt) e-^ t - vx+v2t/2) . (15) 

The solution of eqn. (|13p can be found by using a mechanical analogy. In fact, interpreting x as the time variable and 
ipo(x) as the coordinate of a fictitious particle, eqn. (j 1 3|) formally corresponds to the Newton's equation of motion of 
this particle (of mass M — 1/2), subjected to the force 

F = -£Vo - g^o +1 ■ (16) 

This force can be derived by the potential 



VW) = fv^ + ^C 2 • (17) 



As any motion of a particle in a potential, this is accompanied by the integral of motion that corresponds to its 
mechanical energy 

^= y(^) 2 + ^o) = const. (18) 

Following this mechanical analogy, it is easy to see that a non-trivial motion can take place only if fi < 0, where the 
typical shape of the potential is similar to the one drawn in Fig. Q] 

Notice that a solution is always given by the equilibrium configuration of the potential (|17p . i.e. 



Mx) = \~f) ' (19) 

This solution can be normalized only on a finite volume L, with the dependence of the chemical potential on the 
volume determined by the normalization condition (|14p . i.e. jl = —gL~ a / 2 : given this dependence of the chemical 
potential, the constant solution is simply 

* " it ■ m 

To determine the ground state, we have to compare the GPE energy of this constant solution with the one of a 
localized wavefunction. For this solution, both ipo{ x ) an d ^dxi x ) should vanish when x — > ±oo. This condition fixes 
the constant value of Tt to be zero (notice that this value is not the GPE energy). In this case, the fictitious particle 
takes off from the origin at x = — oo, moving to the right (or, equivalently to the left, since the original equation is 
invariant under ^o(^) — > ~' i Po(x)), until it reaches the inversion point A at the time x = 0. Once the particle arrives 
in A, it inverts its motion and comes back to the origin with a vanishing velocity. It is clear from this analogy that 
A will be the maximum of the bright soliton solution. 

What we said, however, is not the end of the story. In fact, the kind of motion we have just described occurs for any 
potential with the shape shown in Fig. [TJ But we are looking for that particular motion that satisfies the additional 



FIG. 2: Ground state wavefunctions of the GNLSE (|13[) for a — 1,2, 2.5, 3 (dimensionless units are used with g — 1). 



constraint (|14| and this condition can be fulfilled only for a particular shape of the potential, i.e. for a particular 
combination of the parameters fl and g: it is as if the solution is looking for its proper potential. 
In the following, it is convenient to introduce the quantities 



a2s _Mo±2) >0 



2.9 



4g 

a + 2 



7 : 



(21) 



Using the integral of motion Ti., for generic values of fl, g and a (with fi < 0, g > and a > 0) the solution is given 
by a quadrature 



dq 



ijj (x) 



r/r 



(22) 



where ^4 = a 7 is the inversion point reached by the particle at the "time" x — 0. Using the exact expression of the 
integral of the left hand side 



dq 



1 



one gets 



\Ja 2 q 2 - q a + 2 



r: 



a+y/a 2 -4>8(i) 
a-y/a?- iP%(t) 



A 



cosh 7 (%^2pLx) 



(23) 



(24) 



It remains now to impose the normalization eqn. (|T4|) to the solution (|24|) : this fixes the shape of the potential, i.e. 
the relation between fl and g 



(-A)" 



a + 2 



2/Q aT(2/a + l/2) 



27rr(2/a) 



(25) 



When a ^ 4, we can use this equation to express fl as a function of g and, in particular, to write the normalization 
A as 



.4 



/ 2g r( 7 + 1/2) 

tt7(7 + i) r( 7 ) 



(26) 



In Fig. [2] we plot the soliton solution (|24| for different values of a at g = 1. 

However, when a = 4, eqn. (|25p leaves /} undetermined: this means that the corresponding wavefunction 



cosh(2 v / — 2/ix) 



1/2 



(27) 



FIG. 3: Wavefunction (|27)l for N = 3 and g = g* plotted for p, = — 1 (solid line) and fl = —5 (dot-dashed line). Inset: 
corresponding potential (JTTJ) for p, = — 1 (solid line) and /x = — 5 (dot-dashed line). 



is solution of the nonlinear Schrodinger equation (|13|) for every fx. In this case, however, only a particular value of g, 
given by 

9* = , (28) 

guarantees its correct normalization (|14p . Expressed in more physical terms, the attractive 3-body interaction has 
the peculiarity that one can arbitrarily vary the chemical potential provided that the coupling constant be fine-tuned 
to the critical value g*: increasing or decreasing (in modulus) the chemical potential simply results, in this case, in 
shrinking or enlarging the shape of the soliton. This is shown in Fig. [3] where the wavefunction (j2~T)) is plotted for two 
different values of fl: in the inset we plot the corresponding potential (|17p . showing a larger (smaller) inversion point 
corresponding to the smaller (larger) width. 

The fact that fl is undetermined and the soliton can change arbitrarily its shape does not imply that the GPE 
energy is undetermined: in fact, the explicit computation of the next subsection shows that, in this case, the energy 
does not depend on the value of fl. Then, for TV = 3 and g = g*, an infinite degeneracy parametrized by the chemical 
potential fl < occurs. 

In Fig. IHa)-(b) we plot the chemical potential fl for two different values of g, one smaller than g* and the other 
larger: it is seen that for g < g* (g > g*), then fx — > (fx — * — oo) for a — > 4~ while fl — > oo (fx — > 0) for a — > 4 + . The 
singular nature of the 3-body interaction can then be recovered by studying the limit a — > 4 of the formulas (1241) , 
([23]) . (|26[) given above: for a — ► 4T , if g = g* the normalization A goes to 1, while if g < g*, A — > and fl — ► (i.e. we 
have a non-localized solution) whereas if g > <?*, both A and fl diverge, i.e. the wavefunction collapses to the origin. 
It is worth to mention that a singular behavior of the nonlinear Schrodinger equation relative to the value a = 4 has 
been also observed in the dynamical blowing up of the moving wave-packets of this equation: the interested reader 
is referred to the mathematical literature for a detailed discussion of this issue [371 ]. In the present application, this 
instability means that the local 3-body local attractive interactions cannot sustain a bound state unless there is a fine 
tuning of the interaction. In the next Section we will show how an external trap can help to stabilize the bound state 
for a generic value of the coupling. 

To understand better the behaviour of the solution ipo(x) as a function of a, let's define the width a a as cr 2 = 
/ dxx 2 ^{x). One gets 

2 _ £(7+1/2) 7(7 + 1) 
~ 7rr( 7 ) ' 2A« ' a ' [ > 

where X 7 = / dX X 2 / cosh 27 (X). One finds a\{g) = n 2 /3g w 3.28/ g and a 2 (g) = (tt 2 - 6)/(12 5 ) 1 / 3 w 1.69/ g 1 / 3 . For 
large a one has cr 2 — > g/2, while, of course, cr 2 — *■ oo for a — > (no localized soliton without interaction). For g < <?*, 
from eqn. (|29|) one sees that for a — > 4 _ , a Q — > oo, while for a — * 4 + , (Tq, — > 0. In Fig. [5] we plot a 2 for g = 1 < g* 
from eqn. (j2"5|) and, for completeness, also the widths for some value of a obtained from the numerical GNLSE. A 
divergence is observed for a — > 4~, corresponding to the 3-body attraction: the bright soliton becomes larger and 
larger getting close to a = 4, while for a slightly larger than 4 the soliton becomes extremely narrow. This means 
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FIG. 4: (a)-(c) Chemical potential and energy (in dimensionless units) of the bright soliton solution (|24[1 for g = 3 < g* (solid 
line) and g — 4 > g* (dashed line) for a < 4 - (b)-(d) chemical potential and energy for g = 3 (solid line) and g = 4 (dashed 
line) for a > 4. 




FIG. 5: versus the nonlinearity degree a. Solid line: 
ground state of the GNLSE (Q3). The value g = 1< g* 



eqn. (|29[) : open circles: (x ) from the numerical determination of the 
and dimensionless units are used. 



that there is a collapse of the solution (pM)) going to a = 4 from large values of a. At variance, for g > g* , then for 
a — * 4~, er Q — * 0, while for a — ► 4 + , er a — > oo. It should be stressed that, for a > 4, although (f24|) is a solution of the 
GNLSE (|13p . it is no longer its ground state: the divergence of a a for a — > 4~ is signaling the disappearance of the 
bound state due to the 3-body interaction. 



A. Ground state energy 



Using the bright soliton solution we can now estimate the energy per particle. Going back to the physical 
dimensions of all quantities and normalizing now ?/>o to N't, for a ^ 4 the chemical potential is given by 

M=-~ IVm ( 15~7~~T TT ^l ' ( 3 °) 



2m/ 7 2a/(4 " a) Vfi 2 7(7+1) 

where 

V^r( 7 ) 
/7 r( 7 + i/2) 
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FIG. 6: Ground state energy vs. a. Solid line (dashed line): eqn. (|31[l for a < 4 (a > 4); open circles: energy of the numerical 
ground state of the GNLSE (I13[) . Dimensionless units (with At = 1) and g — 1 < g* , as well as different scales of the energy 
for a < 4 (left part) and a > 4 (right), are used. 



The energy per particle is then obtained from the GNLSE energy functional 

2c 



Egp 



Using eqn. (|24l) we obtain 



where 



dx tpQ (x) 



Egp 



h 2 d 2 
2m dx 2 



a - 



At 



h 2 (2mG 
2m 



r, 2 



Mx) r 



£{a) 



ip (x) 



£{a) = 



7(7- 



7 2 / 7 2 



a 



:7(7 + l)/ 7 / 7 +i 



For iV = 2, it is /z 



2 Af^/8h 2 and the previous energy © is recovered. From eqn. 



(31) 



(32) 



(33) 



it follows that in order 



to maintain finite the energy per particle for large At one has to keep G fixed. By a numerical determination of 
the ground state of the GNLSE, we have verified that (|24|) indeed coincides with the ground state for a < 4 both for 
g < g* and g > g*. In Fig. [S]we compare for g — 1 < g* the ground state energy per particle from eqn. (|32p with 
the ground state energy obtained for some values of a obtained from the numerical GNLSE. For g > g* and a < 4, a 
similar agreement is obtained. 

For a = 4, as discussed in the previous section, the chemical potential is undetermined. However, a direct sub- 
stitution of (J57J) in (J3TJ) reveals that E GP = for g = g* . Since (J27J) is a solution of the GNLSE (JT3J) for arbitrary 
fi < 0, and then with arbitrary width, we conclude that an infinite degeneracy - parametrized by a negative chemical 
potential - occurs. 

Using the energy (|3"2"|) we can also estimate the energy of the constant solution in the finite interval [-L/2, L/2}: it 
is E cons t/MT = — 2cp a l 2 j(u + 2), where p = Mt/L is the density. To compare this energy with the previous result 
(|32p for the bright soliton solution we have to choose how to perform the thermodynamic limit: if we choose to keep 

fixed the quantity G — cM^ 2 , with large but finite value of At, sending L to infinite the energy of the constant 
solution vanishes. This means that for a > 4 the constant solution is the ground state of the system; a — 4 and 
g = g* is the case in which both the solutions have zero energy. 



IV. EFFECT OF AN HARMONIC TRAP 



In three dimensions, the presence of an external harmonic trap (with frequency to) can help to stabilize the attractive 
two-body interaction: the critical value of the particle number that induces the collapse, obtained from the GPE, 
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UNSTABLE 




FIG. 7: Stability region according to eqn. (|38fl : D = 1, 2, 3 corresponds respectively to a = 4, 2, 4/3, i.e. iV = 3, 2, 5/3. 



is given by 0.57 N't | a | /a osc [23j]. This critical value can be estimated by a variational method [l6|,[24]]: a gaussian 
trial wavefunction, with the width variational parameter, is introduced and the corresponding energy computed [25j . 
Without external trap (u> — 0), the energy does not have a minimum. However, with lo ^ 0, two situations are 
possible: for a number of particle Nt < N T cr ^ a metastable minimum appears, while for At > N T cr ^ there is no 
minimum. 

In this Section we consider the corresponding problem in one dimension with a 3-body interaction and we show 
that there is a critical value c* of the interaction, such that for c < c* , the bound state is stable. To this aim, we use 
the variational wavefunction 

^ v {x) = C exp(~x 2 /a 2 ) , (34) 

normalized to Nt- The energy to be minimized is obtained by inserting the variational wavefunction (|34j) in the 
generalized GPE functional 



E = / dxtp*(x) 



— — — x 1p(x) +— UJ*X 

2m dx 2 a + 2 1 v w 1 2 



1>(x) . (35) 



To better illustrate the peculiarities of the one-dimensional case, it is useful to perform the analysis also in higher 
dimensions. The energy in D = 1, 2, 3 is given by 

E h 2 „ N$ muj 2 a 2 , , 

where 

a + 2 \a + 2 J \tt J 

Let us consider initially the homogeneous case: for w = the energy (|37j) has a minimum only when D a < 4. For 
Z? = 1, the critical value corresponds to a = 4: this is in agreement with the result of the previous Section, which is 
now obtained by a variational approach. The critical condition 

Da = 4 (38) 

is plotted in FigJT] this figure shows that the higher TV-body interactions tend to be more unstable in higher dimensions. 

When the harmonic trap is present (uj ^ 0), there is still a minimum when the values a is pushed to the range 
of value a < 4/D. When Da > 4, we can identify the critical value c* as indicated in [lij for the D — 3 and 
a = 2 case: since E — > — oo for a — > and E — > oo for a — * oo, the critical value is obtained by the conditions 
dE/da = d 2 E/da 2 = 0. In this way we arrive to the result 



c 



7Vt " I 4mo,2 J {m(Da + 4)J 2a [ tt J \2J ' [ ' 
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FIG. 8: Critical value g* vs. a for a > 4. 
line) and ui = 4 (solid line). 



with D = 1. Dimensionless units (with g = At = 1) are used, with u = 1 (dashed 



The instability curve (|39[) depends on D: with dimensionless units (and At = 1), in one dimension as a — > oo the 
critical value g* goes to zero for u> > n and to infinity otherwise; in two dimensions the behaviour is similar except 
that g* — ► 7r/e when u; = 7r; while in three dimensions critical value goes to infinity for n > lu and to zero otherwise. 
A plot of the g* in D = 1 for a > 4 is presented in Fig. [5] 

Let's now examine the critical point D a = 4, which in £) = 1 corresponds to the 3-body interactions. Minimizing 
the energy ([36]) with respect to a, one finds that there is the critical value 



*AC 



•a/2 



2mf a . 



D 



(40) 



For c < c* there is a minimum and the system is stable, while for c* the energy does not have ever a minimum: 
hence, irrespectively of how large u) may be, the system always collapses. For uj = (no trap) and c < c* (c > c*), 
the minimum value of the energy is then obtained for a = (<r = oo). Notice that the variational approach gives 
for D = 1 and a = 4, the critical value (in dimensionless units) 5* = 3v / 3tt/4 rj 4.08, in good agreement with the 
analytical value g* = 3ir 2 /8 » 3.70. 



CONCLUSIONS 



In this paper we have analyzed the one-dimensional Bose gases with A-body local attractive interactions: by using a 
mean-field approach, we found that A = 3 (i.e., a — 4) is a critical point, and that the localized solution is possible 
only for a critical value c* of the interaction strength. For this critical value, an infinite degeneracy occurs: this 
degeneracy is parametrized by the chemical potential (i.e., eigenfunctions with the same negative fi has the same 
energy). For a < 4, the bright soliton coincides with the ground state of the Gross-Pitaevskii equation: when c < c* 
(c > c*), then a — > 4~ gives a width going to diverge (vanish). We have also studied how an harmonic trap can make 
stable this bound state, pointing out that for A = 3 there exist a critical value of the 3-body interaction strength, 
and below such critical value the localized state occurs. Above this critical value, the collapse is not prevented even 
for very large trap frequency. A brief discussion of the role played by the dimension for A-body local attractive 
interactions has been also presented, showing that higher body interactions are more unstable in higher dimensions. 
We also mention that two-body nonlocal attractive interactions has been studied, showing different ranges of stability 
with respect to the local ones [26| : we could then expect that for 3- and A-body interactions this effect could become 
even more relevant. 

Several proposals have recently addressed the issue of inducing and controlling three-body terms. In [20] it has been 
proposed to use cold polar molecules driven by microwave fields to obtain strong three-body interactions, controllable 
in a separate way from the two-body interactions, which in turn can be switched off [20]. Three-body interactions can 
be effectively induced in mixtures of bosonic particles and molecules: in [l8[ the ground state of rotating Bose gases 
close to a Feshbach resonance has been studied, showing that for suitable parameters they are fractional quantum 
Hall states, whose excitations obey non-abelian exchange statistics. In [19| it was shown that a system of atoms and 
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molecules in a one-dimensional lattice can be effectively modeled by a three-body local (i.e., contact) interaction, 
characterized by a strength U and in the limit U — * oo (without a two-body interaction) the ground state properties 
were investigated by a Pfafhan-like ansatz. The strength U of the three-body interaction can be made also negative 
by using Feshbach resonances. 

One of the main reasons of interest of these proposals relies on the fact that exotic quantum phases, such as 
topological phases, appear to be ground states of Hamiltonian with three or more body interaction terms, an example 
being the fractional quantum Hall states described by the Pfaffian wavefunctions [38[ . The excitations of Pfaffian states 
are non-abelian anyons, on which schemes of fault-tolerant topological quantum computation are based [39j |. We think 
that, in perspective, the possibility to induce and tune effective 3-body interactions could become an important tool 
to control the nonlinear dynamical properties of localized wave-packets and to induce new exotic strongly correlated 
phases in ultracold atoms. 
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